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Abstract. The addition of tunnel barriers to open chaotic systems, as well as 
representing more general physical systems, leads to much richer semiclassical 
dynamics. In particular, we present here a complete semiclassical treatment 
for these systems, in the regime where Ehrenfest time effects are negligible and 
for times shorter than the Heisenberg time. To start we explore the trajectory 
structures which contribute to the survival probability, and find results that are 
also in agreement with random matrix theory. Then we progress to the treatment 
of the probability current density and are able to show, using recursion relation 
arguments, that the continuity equation connecting the current density to the 
survival probability is satisfied to all orders in the semiclassical approximation. 
Following on, we also consider a correlation function of the scattering matrix, for 
which we have to treat a new set of possible trajectory diagrams. By simplifying 
the contributions of these diagrams, we show that the results obtained here 
are consistent with known properties of the scattering matrix. The correlation 
function can be trivially connected to the ac and dc conductances, quantities of 
particular interest for which finally we present a semiclassical expansion. 



PACS numbers: 03.65.Sq, 05.45.Mt 

1. Introduction 

Quantum systems that are chaotic in the classical limit exhibit universal behaviour 
that can be well modelled by random matrix theory (RMT) which involves treating 
the Hamiltonian as a matrix with random elements [1, 2, 3]. Another approach is to 
use semiclassical methods to obtain approximations in terms of classical trajectories, 
which are valid in the semiclassical limit ?i — > 0. These methods, as well as providing 
an explanation of the observed universal behaviour in terms of classical correlations, 
can also be used to explore system-specific quantum properties [4, 5]. 

The semiclassical trajectory based techniques which we use here were first 
developed to treat the spectral form factor K{t), which for quantum chaotic systems 
has a universal form depending only on the symmetries of the system and which can be 
obtained from RMT. The form factor can be written semiclassically using Gutzwiller's 
trace formula [6, 7] as a double sum over periodic orbits. By pairing periodic orbits 
with themselves (or their time reverse), known as the 'diagonal' approximation, it was 
shown [8] , using the sum rule of [9] , that this recreated the leading order term in the 
small r expansion of the form factor. 

Contributions beyond the diagonal approximation come from pairs of correlated 
periodic orbits whose action difference is small on the scale of h. The first such pair 
was found in [10] for a system with uniformly hyperbolic dynamics, and is depicted 
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Figure 1. The type of periodic orbit pair that gives the first ofl'-diagonal 
contribution to the spectral form factor for systems with time reversal symmetry. 




Figure 2. The types of periodic orbit pairs that give the first ofE-diagonal 
contributions to the form factor for systems without time reversal symmetry. 
The encounter regions are indicated by the rectangles. 

schematically in figure 1. It consists of an orbit with a small angle self crossing and 
a partner that follows almost the same trajectory, but which avoids crossing and 
completes the trajectory back to the crossing in the opposite direction. Such pairs 
can therefore only exist in systems with time reversal symmetry, and were shown to 
give the first off-diagonal correction to the form factor, which agrees with the second 
order term of the orthogonal random matrix result [10, 11]. 

Not long after these ideas were reformulated in terms of phase space coordinates 
instead of crossing angles [12, 13, 14], the orbit pairs responsible for the next order 
correction were identified [15, 16], and their contribution shown to agree with the 
next term in the RMT result. Of these orbit pairs, those that are possible for systems 
without time reversal symmetry are depicted in figure 2. 

These ideas and calculations were further extended in [17, 18] to cover orbits with 
an arbitrary number of encounters each involving an arbitrary number of stretches. 
A self-encounter that involves I stretches of the trajectory is called an 'Z-encounter', 
and the encounter stretches are separated by long trajectory stretches called 'links'. 
By using the hyperbolicity and long time ergodicity of the chaotic dynamics, as well 
as considering the number of different possible configurations of orbit pairs, they were 
able to generate all terms of the small r RMT expansion for the unitary, orthogonal 
and symplectic symmetry classes [17, 18, 19]. Since then, these methods have been 
successfully applied to show agreement with RMT, for t < 1, for the transition between 
the unitary and orthogonal symmetry classes, parametric correlations, open systems 
and combinations of all of these [20, 21, 22, 23]. 

Recently though, exciting progress in the regime t > 1 has come from these types 
of methods by considering a generating function of the correlation function using 
spectral determinants [24, 25]. This, along with the use of resummation [26] which 
re-expresses the sum over long periodic orbits in terms of shorter ones, allowed the 
correlation function to be expressed in terms of a sum over four sets of periodic orbits 
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(or pseudo or composite orbits) and then the full form factor for all r to be recreated. 

On a different front, the application of these methods to quantum transport 
follows a similar history as for periodic orbit correlations. One quantity of particular 
interest is the conductance [27], which is given semiclassically by trajectories which 
start in one lead and travel to another [28, 29, 30]. The diagonal terms were evaluated 
first [31, 32], before the first off-diagonal contributions were identified [33]. These are 
related to the periodic orbit pairs in figure 1 and can be formed by cutting the link 
on the left and moving the cut ends to the leads. Identifying the connection between 
the possible periodic orbit structures and the open trajectory structures, and building 
on their work on spectral statistics, the semiclassical expansion for the conductance 
was then calculated to all orders [34]. The same treatment was successfully applied 
to the shot noise [35], conductance fluctuations and other correlation functions [36], 
giving results that agreed with RMT (where RMT predictions existed) and a complete 
semiclassical treatment of quantum transport in that regime. These methods have 
then been applied to include spin interactions [37], to treat the time delay [38], and 
to provide the leading order contributions to higher order correlation functions [39] . 

However, the power of semiclassical methods is not just restricted to recreating 
RMT results, since they can also be applied to regimes where RMT no longer holds. 
One example is the periodic orbit encounters treated in [38, 40] which semiclassically 
recreate the oscillatory terms in the time delay, while a large area of interest is 
in Ehrenfest time effects. For the conductance, the Ehrenfest time dependence 
of weak localisation [41, 42] and coherent backscattering [43, 44] has been treated 
semiclassically. Beyond this, the Ehrenfest time dependence has also been found for 
the leading order of the shot noise [45], conductance fluctuations [46] and a third 
order correlation function [47]. This work has shown that we need to treat additional 
trajectory diagrams that only play a role when the Ehrenfest time is important. 
These include the coherent backscattering contribution, which has an Ehrenfest time 
dependence and can no longer be treated as part of the diagonal approximation, as 
well as the periodic orbit encounters which appear for the conductance fluctuations 
[46, 48]. 

The semiclassical treatment for the conductance and related quantities involves 

trajectories that start and end in the lead, so (apart from coherent backscattering) 
there can be no encounters at the end of the trajectories as the trajectory must escape 
and can no longer return to a nearby point. However, when tunnel barriers were 
included in [49], because the semiclassical trajectories can be reflected when they try to 
escape, encounters can now occur at the leads and additional diagrams exist. Indeed, 
the 'failed' coherent backscattering and other extra diagrams are necessary to preserve 
the unitarity of the evolution [49]. That work was concerned with the leading order 
corrections when Ehrenfest time effects are important, but leaving behind this regime 
similar types of diagrams were also developed to treat the survival probability [50]. As 
the semiclassical approximation for the survival probability involves trajectories that 
start and end inside the system, encounters can occur near the start or the end of the 
trajectory, leading to the 'one-leg-loop' diagrams of [50]. These were generalised in [51] 
where the relationship between these new diagrams and closed periodic orbit structures 
was also explored. These connections were further explored in the context of the 
semiclassical continuity equation in [52] where combinatorial arguments were used to 
show that the continuity equation is satisfled in the semiclassical approximation. This 
is the basis of the current work, and the goal of this article is to combine semiclassical 
trajectory based expansions [10, 17, 18, 33, 34, 36], with the treatment of tunnel 
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barriers [49] and including the extra diagrams and their generahsations [50, 51], to 

provide a complete scimiclassical treatment for chaotic systems with tunnel barriers, 
at least in the regime where Ehrenfest time effects are negligible and for times shorter 
than the Heisenberg time. 



To do this we first examine the changes [49] that adding tunnel barriers brings 
to the system in section 2. We then re-examine the semiclassical continuity equation 
[52], which allows us to explore all the different types of diagrams that contribute, and 
concentrate first on the survival probability in section 3. To simplify the semiclassical 
treatment we shift to the Fourier domain in section 4, and then use recursion relation 
arguments that allow us to combine sums of semiclassical contributions in an elegant 
way. The probability current density, which we treat in section 5, is related to 
the survival probability via the continuity equation, and after again simplifying the 
contributions, we show that the continuity equation is satisfied to all orders with the 
presence of tunnel barriers. In section 6, we consider a correlation function of the 
scattering matrix, for which we have to treat a new set of possible trajectories. These 
have never arisen before in previous contexts but do arise here precisely because of the 
tunnel barriers. Again we can simplify the sum of their contributions using recursion 
relations and we verify that our final results are sensible with a number of consistency 
checks. Finally, we use this to provide an expansion for the conductance in section 7 
before presenting our conclusions in section 8. 

2. Life under a tunnelling regime 

We are considering a chaotic system linked to a lead carrying M scattering channels. 
With perfect coupling to the leads, the survival probability of a long (ergodic) stretch, 
like a link of time t, decays exponentially, e"'**, where /i = Af/Tn is the classical escape 
rate and Th is the Heisenberg time. When we consider an encounter which involves I 
encounter stretches that last tenc each, because the stretches are close and correlated, 
if one escapes they all do, so the survival probability depends only on the time of a 
single crossing of the encounter, e~'^*""' , increasing slightly the survival probability of 
the trajectory as a whole [34] . 

When we include tunnel barriers, the situation changes somewhat [49]. To be 
precise we add a thin potential wall at the end of the lead so that any incoming (or 
outgoing) particle is separated into a transmitted and a reflected part. This simulates 
imperfect coupling to the leads, a situation which often occurs in real physical systems 
from quantum dots to microwave billiards, and so makes the theoretical semiclassical 
treatment of much wider experimental relevance. In particular we take the limit of tall 
and narrow barriers, so that hitting the tmmcls barriers can be treated as a stochastic 
event, where a trajectory has probability pm to pass through on hitting channel m 
(and probability 1 — pm to be reflected). The survival probability of a single classical 
path is then 



As we consider encounters involving I encounter stretches, we also need to know their 
joint survival probability. The probability for all I stretches to survive upon hitting 
channel m is {l—PmY, so the probability to not survive is simply 1 minus this quantity. 



M 




(1) 



m—1 
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leading to an escape rate [49] 

^ 1 - (1 -pmY 

«=/^Z. M • 

m— 1 

With these changes to the survival probability of trajectories, we are now ready to 
treat the semiclassical approximation to the survival probability of an initially trapped 
wavepacket. 



3. Survival probability 

The survival probability of an initially trapped wavepacket is given by 

p{t) - / dr^(r,t)7/;*(r,<), (3) 

J A 

where A is the volume of the corresponding closed system and tp{r, t) is the solution 
of the time dependent Schrodinger equation. As the wavepacket is initially trapped, 
we have p(0) = 1. The survival probability is linked to an integrated current density 
via the continuity equation 

d f 

-pit) + J{t) = 0, J{t) = JJ{r, t) ■ h^x = 0, (4) 

where S is the cross-section of the opening, with n^, is the vector normal to this section 
at the point x in S, and 

j(r,i) = ^[V'*(r,t)VV(r,t) - V(r,i)Vr(r,t)], (5) 

is the probability current density. 

The semiclassical approximation to the survival probability p{t) was treated in 
[50, 51], and we use those results as the starting point for our inclusion of tunnel 
barriers. The diagonal approximation, which only involves a single trajectory that 
starts and ends inside the system, leads to the simple result 

= (6) 

To move beyond the diagonal term we consider trajectories that have close self- 
encounters, working along the lines of [34, 36]. As the survival probability was 
considered in detail in [51], wc briefly point out the main results here as well as 
the necessary modifications due to the presence of the tunnel barriers. Trajectories 
are labelled by a vector v, whose elements vi list the number of i-encounters along the 
trajectory, the total number of which is V = ^ f;. The number of links of the related 
closed periodic orbit is L = ^ Ivi from which we can generate the open trajectories 
by cutting each of those links, meaning that the number of trajectory structures N{v) 
corresponding to the vector v is closely related to the number of closed periodic orbit 
structures. The semiclassical contribution is separated into these three cases: 

A where the start and end points are outside of the encounters, 

B where either the start or the end point is inside an encounter, and 

C where both the start and end point are inside (different) encounters. 

These three cases are illustrated in figure 3 for a trajectory structure with two 
2-encounters. 
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Figure 3. Examples of the three different types of trajectory structure that 
contribute to the survival probability. In (a) both ends are outside of the 
encounters, in (b) the start point is inside the first encounter while in (c) the 
end point is also inside an encounter. 



3.1. Case A 

To generate the trajectory structures for this case, we simply cut one hnk of the 
corresponding closed periodic orbit. This forms two links in the trajectory structure, 
giving a total L + 1 links. For example, if we cut any of the links of the periodic orbit 
in figure 2a, we arrive at the structure in figure 3a. From a trajectory 7, we can create 
a partner trajectory 7' by reconnecting the encounter stretches, leading to an action 
difference of AS* « su, where the vectors s and u contain the appropriate differences, 
along the stable and unstable manifold respectively, of the original encounter stretches. 
The semiclassical contribution of all structures labelled with a common vector v can 
be written as 

p„,Ait)^N{v) y"dsdMw;„,A(s,w,t)e-'^*-Pe-^«", (7) 

where WT,,A{s,u,t) is the weight of such encounters. The exponential term e~^*'=''p 
is the average survival probability of the trajectories and requires the corrections 
described in section 2. We label the V encounters by a, which each involve la encounter 
stretches that last t"^^. We further label the L + 1 links by i, which each last i,;, and 
then the exposure time is given by 

L + l V V 

i— 1 a— 1 a— 1 

To simplify the calculation, we rewrite (7) as 

P„Mt)^N{v) y"dsdMz„^A(s,M,i)e-^^*e^"", (9) 

using an augmented weight, z.v ji^{s,u,t), which includes the term from the survival 
probability of the encounters coming from the right hand side of (8). 

z„,A(s,M,i) =«;.,A(s,M,i)e^=(^i'°-'''=)*o°„c (iq) 

jt - Eg U^ncf Ua (1 + (Ml^c - Mi JCc) 

where the weight comes from [34] and we have expanded the exponent to first order. 
From here the semiclassical contribution can be easily found as it only comes from 
those terms where the encounter times cancel exactly [17, 18]. The integrals over s 
and u in (9) provide a factor of (27r?i)*^'''~^^^^~^-' where / is the number of degrees of 
freedom of the system. This factor can be combined with the phase space volume 
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appearing in (10) and rewritten in terms of the Heisenberg time as Th = Q/{2'Kti)-l ^, 
while the number of structures N{v) can be found in [19]. 

3.2. Case B 

Starting from the structures for case A, we can also shrink either the link at the start 
or the end of the trajectory so that it moves into an encounter. For example, if we 
shrink the first link of the structure in figure 3a, we arrive at the structure in figure 3b. 
This case corresponds to the 'one-leg-loops' (111) of [50], and we write this contribution 
as 

PvAt) = ^i^) J du z^,b{s, u, t)e-i'''e^^^, (11) 

where again, once we have the expression for the augmented weight, we can find 
the semiclassical contribution easily. Here we have L links in total and an integral 
over the position of the encounter a' relative to the start or end point. However, 
this integral can essentially be replaced by the factor of f^^^ following a change of 
variables [51]. When we generate the trajectory structures by cutting the links of the 
corresponding closed periodic orbit we get the encounter a' at the start la' times, and 
an equal number of times at the end. Upon dividing by an overcounting factor L, the 
augmented weight can be simplified to 

, . ,^ ^ 2 (E„ /.Cc) {t - Eg lo.t".nf~' n« (1 + (Ml^g - W J^e"nc) 
^Ms.U^t) X!f|--n„*?nc ' ^''^ 

and we find the contribution as before. 

3.3. Case C 

The final possibility occurs if we additionally shrink the remaining link at the start 
or the end so that both the start and end point are inside different encounters. For 
example, when we shrink the remaining end link of the structure in figure 3b, we arrive 
at figure 3c, and this is the Oil case in [51]. This contribution can be written as 

P^,c{t) = N{v) J ds du z^,c{s, u, t)e-^''*ei^^, (13) 

where we now have L — 1 links in total and two integrals over the position of the start 
and end encounters relative to the start and end point. The number of links of the 
closed periodic orbit, divided by the factor L, which connect encounter a to encounter 
P is denoted Afa,0{v) and recorded in a matrix JV(v), whose important elements are 
tabulated in [51]. We include the number of possibilities with the augmented weight, 
which reduces to 



L-2 



(L-2)!fi^-^n.i?„c 
n(l + (/^i^«-WJCnc), (14) 



from which the contribution can be found as before. 
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3.4- Unitary case results 



By simply adding the results for the three cases we can obtain the results for each 
vector and for each symmetry class. For the unitary case, the first off-diagonal 
contributions come from a vector with two 2-encounters (which we will denote by 
(2)^), which gives a contribution of 



-lilt 



rp2 



1 2 



^2 ^3 _|_ (2/^1 



M2j 



2 3 24 

while a vector with a single 3-encounter. denoted (3)^, provides 



-flit 



rp2 



1 2 



M3 ^3 



(15) 



(16) 



2 6 

We can sum over the different possible trajectory structures, including these ones, and 
obtain the following expansion for the survival probability 



pit) 



-Hit 



1 



/XI - 2/i2 + , (2^1 - /i2)^ t 



Ml 



6 

4/X2 



t2 



+ 



2/X3 



24 
4/X4 + Ms 



rp2 



15 

13/ii/i3 - 3/U1/U4 



7/i2M3 



+ 



29/ii/X2 



60 



30 y 



(17) 



45 36 

However, in order to simplify this general result, we can set all of the individual 
tunnelling probabilities to p. The Heisenberg time dependence involves only the value 
of L — F of the vector, and so we further set t = tTh- When we do this for the unitary 
case we obtain 

p2(p -1)M^3 ^ p^ip - 1)M_5 , p^9p - 7)M2_g 



Pit) 



-pMr 



24 



15 



180 



p^{7p-3)M^ p^ip- 



+ 



720 2, 
p^M^ /(401p- 356)M2 



my. 



(18) 



1920 10080 

which reduces to the previous result [51] upon removing the tunnel barrier (by setting 
p = 1). In this form it is easier to see the effect of changing the tunnelling probability, 
especially if we keep the escape rate (or pM) constant. The first off-diagonal term, 
which is due to the interplay between trajectories with two 2-encounters and those 
with a single 3-encounter, actually increases as the tunnelling probability is decreased 
from 1 (to 1/2) before falling back to again. This is due to the fact that the 
survival probability of a encounter with three correlated stretches falls more slowly 
than that of an encounter with two stretches leading to an overall increase in the 
survival probability. With all the tunnelling probabilities equal, we can compare our 
semiclassical result with the random matrix result of [53] and find full agreement. The 
random matrix result is also completely general, with different tunnelling probabilities 
for each channel, but for the comparison with our semiclassical result we need to 
perform the integrals in the random matrix result and this only becomes feasible when 
we set all the tunnelling probabilities equal. The semiclassical result here in (17), of 
course, provides a direct way of obtaining an expansion for the survival probability 
also in the case where the tunnel probabilities differ in each channel. 
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3.5. Orthogonal case results 

For the orthogonal case, the first off-diagonal terms comes from trajectories with a 
single 2-encounter [10, 33], which give a contribution of 



Pi2v{t) = ^ 

The next order terms again come from structures with two 2-encounters 



P(2)2(i) 



-flit 



rp2 



2t 



2 5(2/11 - 112)^3 ^ 5(2/xi - )i2) 



24 



-t" 



and with a single 3-encounter 



P(3)i(i) = 



-Ml* 



rp2 



.2^2 + 6Mi-2M3^3- 



(19) 



(20) 



(21) 



Summing over different structures, wc obtain the following expansion 



P{t) 



-Hit 



2ni - H2 t^ _ 4Ati - S/tz + t^ 5(2/n - /ta)^ t^ 



+ 



2 Th 3 
12/ii — 23/U2 + I8/U3 — 5/i4 



Tn2 



24 



n-i2 



6 T3 

74/tf + 4:1/4 - 119M1M2 + 30/ti/t3 - 15/X2/ X3 
30 

96/xi - 250At2 + 297/^3 - I68/X4 + 37/X5 



T3 



30 



t-t4 



+ .. 



(22) 



We can again set all the tunnelling probabilities to p and simplify the expansion 
for the orthogonal case, obtaining 



p{t) = e 



_ ^-pMt 



1 



p^M 2 p2(2p-l)M 3 /S/M^ p^{5p'^ -2p-l)M 



-T — 



24 



p'*(15j3-4)M2 p2(-37p3_l7p2_5p_3^^^^ 

30 30 
41/M3 p^{61p'^ - 22p - 4)M2 



720 
+ 



60 



p2(337p4 - 152p3 _ 5ip2 _ 22p - 16)M\ 



180 



/(182p - 37)M3 _^ /(2408p3 - 936392 - 231p - 65)M2 



+ 



840 

35^(1888/ - 877p^ - 266p^ - 129p2 - 76p - 60)M\ 7 
630 ^ 



1260 

„2 



+ ... (23) 

which again reduces to the previous result [51] and agrees with the random matrix 
result [54] , if wc transform it following the steps in [23] . Here wc can sec that reducing 
the tunnelling probability (at fixed escape rate) leads to a direct reduction of the first 
off-diagonal term as the enhancement of the survival probability originally due to the 
closeness of the encounter stretches in the 2-encounter becomes damped. Without 
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tunnel barriers, the higher order corrections are all due to the slight enhancement to 

the survival probability that having encounters brings. But as the tunnel probability 
is reduced to (at fixed pM), this advantage, along with the off-diagonal corrections, 
vanishes. 

Using the same semiclassical techniques, we can also treat transport quantities 
like the conductance and the reader interested in those results might skip straight 
to section 7. Before we arrive there though, we show in the next three sections how 
physical properties like continuity arise from semiclassical recursions and how decay 
is related to transport. To do this we first move to the Fourier space and consider the 
survival probability and later its connection to the current density. 



4. Transformed survival probability 

We start with a trapped wavepacket {p{0) = 1) and restrict ourselves to positive times, 
so we examine the (one-sided) inverse Fourier transform of the survival probability p{t) 

P{uj)= / dTp(rrH)e2---, (24) 



where we still have r = </Th. As we saw in [52], the semiclassical contribution 
for this transformed survival probability can be separated into a simple product of 
contributions from the encounters and the links, as is possible for the conductance 
[34]. For example for case A the contribution from structures corresponding to v is 

P„,,.).^(];]'jf%,e-(».-*)'.) ,25, 

enc f 



11 / 'isadua QlZ^ 



enc 



where we have separated the exposure time using the first expression in (8). When we 
evaluate the integrals, the Heisenberg times cancel meaning that we essentially just 
obtain a factor of (Gi — 27ria;)~^ for each link and a factor of —{Gi^ — 2'r:iu}la) for each 
encounter, where we define 

M 

=MiTH= ^ [1-(1-P™)'] , (26) 
in line with (2). The contribution then becomes 

P.,M=m.)(-i)''n^2^^. ,27, 

Likewise, the diagonal term which involves a single link simply gives 

as can be obtained directly from (6). 

For case B we have one link fewer, leaving L in total, and one encounter, a', 

at the start or end of the trajectory pair. This encounter, which occurs 21^' times 
(divided by L), just gives a factor of 1 leading to the simplified contribution 

PvM^) ^ (g;— 2^kJU) ) (Gr-27ria;)^ " ^^^^ 
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For case C the two encounters at the end both give factors of 1 while the remaining 
encounters and Unks give their usual contributions, providing 

[i^ff ~ 27riu;/„)(G;^ - 2niujlp) j (Gi - 2niw)^-^ ' ^ ' 

4.I. Recursion relations 

To proceed we consider how we can re-express the sum over vectors with a common 
value of L — V for case C using recursion relations. As can be seen in (30) it is the 
sizes of the encounters which is important so we replace Afa.fj{v), by J\fk.i{v) which 
is the number of links connecting a fc-encounter to an Z-encounter in all the closed 
periodic orbits structures described by v (and divided by L). With this replacement, 
each vector then gives the following contribution to the survival probability 

P r-) ^ njGi. - 2niujQ 

If^ (Gfe-27ria;fc)(G,-27ria;Z) I (Gi - 27ria;)^-i ' ^ ' 

Based on the recursion relations in [18, 19], the following relation was obtained [52] 

Afk,i{v) = '''^2'_\^^ N{vi'^^-''\ (32) 

where k' = k + I — 1, and Vk' is the (fc')-th component of v. This relation 
represents merging a k- and Z-encounter in v to create a fc'-encounter in 1 
which is correspondingly formed by decreasing the components Vk and Vi by one while 
increasing the component Vk' by one (so that Vk' + 1 = v^,' ' for example). When we 
include the extra factors to match the form of the survival probability, (32) becomes 

{Gk - 2'Kiujk){Gi - 271^1) (Gi - 2^0;)-^-! 

°- (a!-2llo "'°""'"''"'' 

where 

We now sum (31) over all vectors with the same value oi L — V = m, and use 
the recursion relation (33) to re-express the m-th order contribution to the survival 
probability for case C as 

L — V—m , f f 

V k,l ^ " ' 

where v' = dI*^-'^*^]. To form v' from v we combined a k and Z-encounter which 
reduces both L and V by one, but leaves the value oi L — V = m unchanged. The 
sum over v can then effectively be replaced as a sum over v' itself [18, 19]. We later 
identify this dummy sum variable v' with v when we combine this contribution with 
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the contributions from the other cases to obtain (38). We first recall the contribution 
of case A from (27) 

L-V=m 



L 



(Gi - 27ria;) 



N{v,G), 



(36) 



and from case B from (29), where we re-express the sum over a in terms of the 
components of the vector v 

L-V=m 



2lvi 



I 



{Gi - 27riw/) 



N{v,G), 



(37) 



so that we can express the total contribution to the survival probability as 

L-V=m 



L 



(Gi - 27riw) 



E 



2lvi 



k'vk' 



{Gi - 27riw0 ^ (Gfe/ - 27ria;fc' 



N{v,G). 



(38) 

Moreover, we can simplify the double sum in the third term. Ks k' = k + I — 1 and 
/c, Z > 2 for each value of k' we get {k' — 2) copies of that term. The double sum can 
then be simplified as 

k'{k' - 2)<, 



E 



k'v',, 



E 



^ ^ {Gk' - 27ric^fc') {Gk' - 27ricc;fc') ' 



(39) 



We note that the term /c' = 2 in resulting sum is 0, so we can lower the limit of the 
sum accordingly. By identifying k' with I, we can combine this sum with the sum over 
I in (38), which becomes 

L-V=r 



Pm{^)= E 



L 



(Gi - 27riw) 



E 



(Gi - 2mujl) 



N{v,G). (40) 



This result means that we can effectively replace all of the different types of 
contributions by this simplified form. Of course the contribution of individual vectors 
differs, but for the sum over all vectors with a common value L — V = m, which is 
the more useful semiclassical quantity, this is a very helpful simplification. 



5. Transformed current density 

The current density, which is connected to the survival probability through the 
continuity equation, is expressed semiclassically in terms of trajectories that start 
inside the cavity but end in the lead. Without tunnel barriers, as soon as the trajectory 
hits the lead it escapes so, as we saw in [52], we cannot have case C and we only get 
half the contribution for case B as the encounter can no longer occur at the end. With 
tunnel barriers these cases arc again possible, as long as the trajectory is reflected on 
all but the last encounter stretch in the lead. 

We consider the Fourier transform of the integrated current density J{t) 

J{uj) = / dr J(TTH)e2--^ (41) 
Jo 

for which we can again write the semiclassical contribution as a product of 
contributions from the links and encounters, as for the survival probability. Now 
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however, the end of the trajectory must hit the lead and escape. Whichever channel 
m it hits, it only has a probability p„i of escaping, and this leads to a channel factor 
of Gi = 'Ylim=iPm- For the diagonal approximation we then get 

while for trajectory structures described by v for case A we obtain 

JvA^) = ^PvA^)- (43) 

With case B wc have to treat the case where the end of the trajectory is during 
an encounter differently from the case when the start is during an encounter. When 
the encounter occurs at the start, which happens half of the time, we can treat it as 
we did for the survival probability as long as we include the final escape probability 
factor Gi. When an Z-encounter occurs at the end it must be reflected the first {I — 1) 
times so as to return to and create the encounter, while escaping on the last encounter 
stretch. The probability of this happening for a particular channel is then simply 
(1 — Pmy~^Pm, which we must sum over all channels and include as a factor too. If 
we define 

M 

Hi=J2Pmi^-Pmy-\ (44) 
m=l 

then the total contribution can be expressed as 

JvA^} - j;^ (G;^, -2W„0 ) (Gi-27ric.)i " ^^^> 

For case C we know we have an encounter at both ends so, if we say that encounter 
a is at the lead, while /? is at the start of the trajectory, we can simplify the result to 

J r(co) = f V ^^Av)Hi. ] (-l)^-^n.(g^.-27ria;U 

^^'<' - '^^''^^c.m, - 27ric^Z^) j Tu{Gr - 27ria;)^-i ' 

Note that the matrix A/'a,/3(v) and hence the double sum is symmetric under swapping 
a and /?. 

When we combine the different cases, and re-express the contribution from case 
C using the recursion relations in section 4.1, we obtain 



EiL^LTi -r ni) X - 
^ {Gi - 27ria;/) ~ ^ {Gk' - 2mu>k') 



N{v,G) 



(47) 

Again we can simplify the double sum in the third term. For each k' we obtain a copy 
of each Hh for fc = 2, . . . , /c' — 1. The double sum is then 

fe,J ^ ' k'>2 

The k' = 2 term would involve no sum over Hk and is formally zero so we can again 
lower the limit of the sum. We can then combine this sum with the sum over I in (47) 
which becomes 



^ {Gk' - 2TTiwk') {Gk' - 2TTiwk') ' 



L-V=m 



LG^ _ lvi{Gi+Y!k=2Hk + Hi) 



{Gi - 2niuj) ^ {Gi - 2^1(^1) 
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M 



M 



(1 -Pr7 



= G, 



(50) 



k=2 m=l k—1 m=l 

as the sum just involves a geometric progression. The final result for the current 
density is then 

L-V=m 



E 



LGi 



(Gi - 27riw) 



E 



IviGi 



{Gi - 2iTiujl) 



N{v,G) 



(51) 



5.1. Continuity Equation 

We have examined the current density because it is connected to the survival 
probability via the continuity equation (4), which in the Fourier space becomes 

Th J(w) - (27ria;)P(a;) = 1, (52) 

where semiclassically the 1 comes from the diagonal terms, which can easily be checked 
using (28) and (42). As such, to ensure that the continuity equation is satisfied 
semiclassically, we need to show that the off-diagonal terms vanish and that 

Th Jm(a;) - (27ria;)P™(a;) = 0, (53) 

for all m > 0. Combining (40) and (51), we have to evaluate 

L-V=m 



E 



L(Gi - 27ria;) ^ lvi{Gi - 2-^1^1) 



{Gi - 27ria;) 



{Gi - 27ria;/) 



N{v,G), 



which directly reduces to 

L-V=m 

E 



L 



Ivi 



N{v,G) = 0, 



(54) 



(55) 



since Ivi = L. This verifies (53) and that the scmiclassical expansion respects 
the continuity equation. Adding the diagonal terms, we indeed obtain (52) in our 
scmiclassical regime. 



6. Connection to transport 

A large area of interest in semiclassics is the treatment of quantum transport, 
rather than decay, but we have followed this route because, as we saw in [52], the 
current density is connected to a transport quantity F{t) via the continuity equation 
F{t) + dJ{t)/dt = 0. In the Fourier space this continuity equation is 

TyiF{uj) - (27riw) J(w) = 1, (56) 
where the semiclassical approximation to F{lo) is given by 

n^)«^E E i^.i?yeH(^.-^y)e^(*^+V), (57) 

^ a,h 7,7'{a^ti) 

where the sum over a and h is over the channels in the lead and we sum over trajectories 
7 and 7' connecting these channels, which have actions S^, stability amplitudes D-y and 
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times t^. This is a semiclassical approximation to a correlation function of scattering 
matrix elements 

(58) 

which was considered in detail in [38] and is related to the Wigner time delay as 
well as the ac conductance [55] and hence the conductance. We first calculate the 
contributions and then we will later provide an expansion for the conduc;tance. 
The simplest contribution is the diagonal approximation which reduces to 

F'^-^H^^J] J2 iD.l'e'^'^ (59) 

H a,b -Yia^b) 

where we will use the sum rule [33] which has been implicit in our previous calculations 
^ \D,f...^ / dt.e-'^^^... (60) 

7(a^6) 

and modified to represent the survival probability now that we have tunnel barriers. 
We also need to perform the sum over the channels a and b. We remember that we 
have a probability of Pa of tunnelling into the cavity through channel a and likewise 
a probability Ph of leaving through channel b. The sum over channels is then simply 
X^o bPaPb = Gf. Substituting into (59) and performing the integral we obtain 

^""M'=51(g7^' '''' 

where we recall that Gi = /^i/Th- We can consider that we have an additional 
contribution for systems with time reversal symmetry when the start and end channels 
coincide {a = b). Then we can also compare the trajectory 7 with the time reversal 
of its partner 7' giving an additional p'^ for this channel combination. This extra 
possibility corresponds to coherent backscattering (cbs) and must be considered more 
carefully when Ehrcnfcst time effects arc important since we actually have a 2- 
encounter in the lead. It suits our purposes here to include this contribution with 
the other contributions involving a 2-encounter, which we explore in section 6.1, and 
not with the diagonal approximation. 

We can find the contribution of correlated trajectories using the open sum rule 
and an auxiliary weight function as before [34, 36] sphtting the contribution into 
encounters and links as for the survival probability in section 4. For case A, the 
contribution of each vector can be written as 

LGj N{v, G ) 
{Gi - 2nioj) T| 

For case B, we have an Z-encounter at the start or the end so we get a channel factor 
of G\Hi, and an additional factor of 2, giving 



F..A(a;)= ,^ (62) 



^ ^ 2lviG,Hi N{v,G) 

F„,b(w) = -2^— TT-r-j^^pi — ' (63) 

while for case C, with an l-encounter at the start and a fc-encounter at the end, we have 
a channel factor of H}.Hi. To simplify matters we will also sum over all vectors with 
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the same value oi L — V = m and re-express the result using the recursion relations 
in section 4.1 to obtain 

F„,cM = - (G,, - 27ric.fcO Tl ' ^^^^ 

where again k' = k + I — 1. As wc also saw in section 4.1, wc can simphfy the sum 
over k, I and combine the result with the sum over / of the contribution of case B, to 
obtain a total result for cases B and C of 

p ( ^ ^-^"V IviKi N{v,G) 
where we have defined 

I 

ki=Y,HkHi_k+i, (66) 
fe=i 

and where the k = 1 and k = I term effectively come from case B (as Gi = Hi), while 
the rest come from case C. 

However, for transport quantities, where wc can start and end in the same channel, 
there is an additional possibility: Case D, where both the start and end point are inside 
the same encounter (cbs). 



6.1. Case D 

When the start and end channel arc the same, wc have the additional possibility that 
the trajectory can start and end in the same encounter at the lead. The coherent 
backscattering contributions included here which involve a 2-encounter at the lead 
could be included in the cases above, as for systems with time reversal symmetry wc 
can always add a 2-encounter near the lead without affecting the rest of the trajectory 
(though for cases B and C this creates a more complicated encounter at the lead). 
Although they arc included here for practical reason, wc arc particularly interested 
in the possibilities that can only occur with tunnel barriers. These possibilities can 
also occur for systems without time reversal symmetry, marking their difference from 
coherent backscattering, though they are in another sense a generalisation of that case. 
They are also a Oil contribution like case C, but they only involve a single encounter 
at the lead like case B. 

To calculate their contribution, wc remember that for case C wc started with 
the periodic orbit structures and counted all the links that connected two different 
encounters. Obviously, if the link does not connect two different encounters it must 
connect the encounter to itself, and for all structures of type v wc record these numbers 
in a vector N(t;), where the component N;(t)) is the number links connecting an l- 
encounter to itself. We also need to know the channel factor for such an /-encounter in 
the lead. Wc first enter the channel, a, with probability pa, then get reflected {I — 2) 
times before escaping on the last stretch, giving a total factor of /;_!, where 

M 

ii = j2pii-^-pmy-'- (67) 

m=l 
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There are L — 1 links in total, and a more complicated integral over the encounter, 
but the result reduces to 



A,,D(a;)= ij2 



i-i 



27ria;L 



(G( - 2mLol) 



(Gi - 2mLo) 



L-l 



(68) 



Because a link either connects an encounter to itself or a different encounter, we 
have the relation 



N{v) = Y,^iiv) + 



(69) 



k,l 



However, we can, using the relations in [18, 19], break this sum up into contributions 
coming from each I 



Ivi 



'JlN{v)=-Ni{v) + ^Afk,i{v), 



(70) 



where we note that A4,i is a symmetric matrix, and if we sum this result over 
/ we recover (69). We also note that for systems with time reversal symmetry 
J^2{v) = N{v^'^^^) where t;[^^I is the vector formed by removing a 2-encounter from v. 
This relation is the reason why we could also include coherent backscattering involving 
a 2-encounter at the lead as an additional channel factor in cases A, B and C, but 
we instead include this contribution here precisely because of (70) as we already have 
recursion relations for Mk,i- When we substitute the result (70) into (68), and use the 
recursion relation (33) we obtain a contribution of 



E 

L / 



lviIi-i{Gi - 2'Kiuj) 
(Gi - 27riw/) 



N{v,G) 



E 

k,l 



k'v'^,Ii-i{Gk-2mujk) 
{Gk' - 2mu;k') 



Niv',G) 



(71) 



where we recall that k' = k + I — 1 and v' = v^'^'^^'' K When we sum over all vectors 
with a common value of L — V = m we can re-express the second line, perform one of 
the sums and combine the result with the first line to obtain 



-Fm,D(w) 



L-V=m 

E 



E 



2ki^(1 - k)] 



[Gi - 2-k\ujI) 

Using the results which are shown in Appendix B, namely 



N{v,G). 



l-l 



l-l 



- k)Ik = /Gi - Gi, ^4Gi_fc = GiG, - ki, 
fe=i fe=i 
we can simplify (72) to 



-Fm,D(w) 



L-y=m 

E 



E 



Ivi 



GiGi-Ki-2Trhj{lGi-Gi) 



[Gi - 2-k\wI) 



N{v,G). 



(72) 



(73) 



(74) 



Finally, we can combine this result with the other cases in (62) and (65), to obtain 



L-V=m 



LGl 



(Gi - 27ria;) 



-E 



Ivi [GiGi - 27riiv{lGi - GQ] 
(G; - 27ria;0 



N{v,G) 



rp2 



(75) 
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6.2. Consistency 



As a first check of our results for F{u>) in (61) and (75), we consider the case when 
u) = 0, for which 

From the diagonal term in (61), we obtain 



(76) 



F^'^s(o) = 



rp2 ' 



while for the off-diagonal terms we have 

L-V=m 



LGi-^lviGi 



(77) 



(78) 



which is identically as L = '^ilvi. This means that scmiclassically wc have 
Tr [^SS^ = Gi, which is the effective number of open channels, and the result we 
would expect. 

Returning to the continuity equation (56), we can now show that 

TuFmito) - {2TriL0)Jm{uj) 

L-V=m 



E E 



lg,-y: 



IviGijGi - 2Triujl) 
(G; - 2mLul) 



Niv,G) 
Th 



0, 



(79) 



for exactly the same reasons as above and in (54). The result in (79) for all m > 
combined with a direct check of the diagonal terms from (42) and (61), shows that 
the continuity equation (56) is indeed satisfied in the semiclassical approximation. 
As a final check, we know that F{u)) is related to the time delay, as in [38], via 

If we put in the diagonal approximation result from (61) into this equation we obtain 

Gi ^1 ' 

which is exactly the average dwell time of the system with tunnel barriers. For the 
off-diagonal terms from (75) wc obtain 



diag 

'w 



(81) 



Tm,W = 



1 

Ml 



L-V=m 

E 



+ 



Ml 



^+E 

m 

E 



lvi{lGi-Gi) 



L-V=m 



Gi 

LGi-Y,lviGi 



(-i)^aGr 

L Gf 



N{v,G) 



0, 



(82) 



iv=0 



again for the same reasons as above. We therefore find that the average value of our 
semiclassical expression for the time delay docs give the average time delay. Now that 
we have shown that the semiclassical results satisfy unitarity and we have a firm footing 
for the contribution of all the possible cases for chaotic systems with tunnel barriers, 
we can use these results to calculate a semiclassical expansion for the conductance. 
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7. Conductance 

Up until now we have considered a cavity with a single lead containing M channels. For 
the conductance we can imagine splitting this lead into two parts, a left lead containing 
Ml channels, and a right lead containing Mr channels, where Ml + Mr = M. We 
will consider an ac conductance given by 

a— 1 0=1 

where Tba are the elements of the scattering matrix that connect the left lead to the 
right lead, and we have set e = 2nu). The semiclassical approximation is similar to 
that for F{u)), and is given by 

AfL Mb. 

^(^)-^EE E i5.^ye*(^--^v)e^(*^+V), (34) 

^ a=l 6=1 ^-,7'(a^b) 

where the only difference is the change in the channel sum. The factors for each 
channel are as before, but the sum is different, so we define 

AfL A/l 

GL = ^Pa, Hj^,l = ^Pa{l-Pay-\ (85) 

a=l a=l 

and similarly for the right lead. For the diagonal approximation, the channel sum 
therefore gives a factor GlGr and we can write down the result simply as 

C'^He) « j^^y (86) 
(Gi — le) 

where we note that because we start and end in different leads there can be no coherent 
backscattering contribution, and no contribution from case D. For the rest of case A, 
the contribution for each vector follows directly as 

C.Ae) = ^^^N{v,G). (87) 

For case B, if the /-encounter occurs at the start, in the left lead, wc get a factor 
-ffL,iGR, while if it occurs in the right lead we have the factor GL-ffR,;, giving a total 
contribution of 

C„,B(e) = - 2^ _ .^^^ -N{v, G). (88) 

For case C, with an /-encounter at the start and a fc-encounter at the end, we have a 
channel factor of -ffL^fe^R,;, giving a contribution of 

J^kAv)HL,kHn,i \ (-l)^^^n.(G,,-ieU 
= I^E (G,-iefc)(G,-ieO ) (G, - ie)-- " ^''^ 

With these three cases we can find the semiclassical expansion for the conductance 
for both symmetry classes. 

Further, if we sum over all vectors with a common value oi L — V = m for these 
three cases we can simplify the result to 

r ( \ - LKi _ \- IvjKi 



N{v,G), (90) 
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following the reasoning in section 6, and where for the conductance we define 

I 

Ki=Y,HL,kHR,i-k+i, (91) 
fe=i 

and wc note that Gl = ^^l,! a-nd that Ki = GlGr. The extremely simple form 
of the semiclassical result in equation (90) is one of the most important results of 
this paper. This is the reason why we have re-explored the survival probability, with 
its three possible cases and its link to the current density via a continuity equation 
[52], for the situation where we have tunnel barriers. For the conductance we have 
exactly the same three cases, and the recursion relations that we needed to develop 
for the continuity equation can be applied directly to obtain the simple result (90). 
Furthermore, this also simplifies our calculation of the semiclassical expansion of the 
conductance. 



7.1. Unitary case 

For the unitary case, the first off-diagonal contributions come from a vector with two 
2-encounters 

GlGr(G2 - 2ie)^ (gL,2gR + GLgR.2)(g2-2ie) , gL.2gR,2 

= (Gi-i6)5 (GTTk)^ + (Gi-ie)3' (^2) 

and with a single 3-enc(ninter 

GlGh(G:. -3ie) iH^.,(ni + C^LgR.s) 

= (Gl - ier + (GT^^ • 

We can then sum over different possible trajectory structures, including these ones, 
to obtain the expansion 

^ K, Ks _ Ki{Gs-3ie) + K2{G2-2ie) 

(Gl - ie) (Gl - ie)3 (Gi - ie)^ 

Ki{G2-2ie)^ + 8Ks 
(Gi-ie)5 

_ 8Ki{G5 - 5ie) -|- 8K2{G4 - 4ie) -|- 12iV'3(G3 - 3ie) + 16/^4(^2 - 2ie) 

(Gl - ie)6 

12A:i(G3 - 3ie)^ + 24Ki{G2 - 2ie)(G4 - 4ie) + 28/^2 (G2 - 2ie)(G3 - 3ie) 

(Gi-ie)7 

21j^3(g2 - 2ie)2 + 180i^7 , 
^ (Gi-ie)7 ^^^> 

This result has not yet been obtained using RMT, but we can compare to previous 
results if we consider the conductance, which is given by setting e = 0. For the 
conductance, only the first off-diagonal term is known from RMT [56] and this is in 
agreement with the result here. To simplify the result further, we set the tunnelling 
probability in each channel equal to p to obtain 



+ 



G(0): 



(95) 



M 

which agrees with the result in [34] when we set p = 1. Along with the obvious 
reduction of the conductance as the tunnelling probability is reduced, which reflects 
the decreased possibility of entering the cavity in the first place, there is also an 
additional reduction from the higher order terms due to the intricate balance between 
the survival probabilities of all the different sized encounters. 



Semiclassics for chaotic systems with tunnel barriers 21 
7.2. Orthogonal case 

For the orthogonal case, the first off-diagonal contribution comes from a vector with 
a single 2-encounter 

^(2)^ (^) = (Gi - ie)3 + ^T^kp ' 

while the next order terms are 

^(2P(^) = ^T^lip ^731^)4 + (Gi-ie)3' (^^^ 

and 

4G lGr(G3 - 3ie) 4(gL,3GR + Gl^ r.s) 
(Gi - ie)4 + (Gi - ie)3 

When we then sum over different possible trajectory structures, we obtain the following 
expansion 

Ki K2 Ki{G2 - 2ie) - 4K3 



G(3)i(e) = 77T T^^i ^ ?7T JT^a • (98) 



G(e) = 



(Gi-ie) (Gi-ie)2 (Gi - ie)3 

4:Ki{G3 - 3ie) + 5K2{G2 - 2ie) - 2OK4 



+ 



+ 



(Gi - ie)4 

5i<:i(G2 - 2ie)2 - 20^1 (G4 - 4ie) - 24/^2 (G3 - 3ie) 
(Gi - ie)5 

3(i/v:s(;G2 - 2i() - UcS/C-s 
(Gi -ie)5 

60A:i(G2 - 2ie)(G3 - 3ie) - 148j^i(G5 - 5ie) + 41A:2(G2 - 2ie)^ 

(Gi - ie)6 

168ii'2(G4 - 4ie) + 228i^3(G3 - 3ie) + 336K4{G2 - 2ie) + 1348i^6 



(Gi - ie)6 

+ . . . (99) 

This result is also novel and the conductance can again be obtained by setting e = 0. 
For the conductance, we can compare the first off-diagonal term with the random 
matrix result from [56], as well as with the semiclassical result of [49], and find 
agreement. We can look at the simplified case where the tunnelling probability in 
each channel is equal to p and obtain 



G(0) = 



MlMr 
M 



p 3p-2 5p^ -4 I V - + 8p - 12 
21^4 - 4p3 - 36p2 + 84p - 64 



(100) 



43p5 - 18^^^ + 88p3 - 528p2 + 896p - 480 



which also reduces to the resiilt in [34] when wc; remove the; tunnel barriers by setting 
p = 1. Considering the conductance as the tunnelling probability is reduced (at 
fixed pM), we can see that the dip in the conductance from the first off-diagonal 
term also reduces. This dip compensates the enhancement to the reflectance due to 
coherent backscattering and, as noted in [49], as the tunnelling probability decreases 
these trajectories trying to return to their starting lead are increasingly likely to be 
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reflected back into the system and contribute to the conductance instead. At higher 
orders the mechanisms are more compUcatcd, and with new coherent backscattering 
possibilities (case D) can lead to changes in sign, but the overall reduction of the terms 
(as p ^ 0) similarly represents that as trajectories become less likely to leave each 
time they hit a channel, they become more likely to leave through cither lead. 

It is worth remarking here that in the small p limit, in particular where pM < 1 
the average trajectory spends longer in the cavity than the Heisenberg time, as can be 
seen from (81). The semiclassical treatment presented here is then no longer complete 
in that regime, but by the same token, adding tunnel barriers leads to an interesting 
example of where Heisenberg time eff'ects become important in open systems. Such 
systems could therefore be useful for exploring classical trajectory correlations beyond 
the Heisenberg time. 

8. Conclusions 

The main focus of this paper has been on the combinatorial relations that connect the 
different possible trajectory structures and their contributions. These build on the 
connections in [52], where without tunnel barriers we move in a strict hierarchy from 
a scattering matrix correlation function involving only case A, through the current 
density where we add case B, and finally to the survival probability where all three 
cases are allowed (A, B and C). Upon the addition of tunnel barriers, this hierarchy 
is equalised and all three cases contribute for all quantities. Moreover, when we start 
and end in the same channel, a fourth and new type of contribution appears: type 
D which is a generalisation of coherent backscattering. This type of structure was 
previously ruled out, but from (69) for example we can see how it naturally completes 
the set. 

Along with the combinatorics for the trajectory structures, which can all be 
generated from the related closed periodic orbit structures by cutting appropriate 
links, the tunnel barriers add probabilistic combinatorics from the channel factors 
and escape probabilities [49]. We have shown how these all combine and simplify 
mathematically, in non-obvious ways, so that through the semiclassical approximation 
they give results that respect required physical properties such as contimiity. Simple 
physics is mirrored by the complicated combinatorial structure, here made more 
complicated by the tunnel barriers. The semiclassical methods [34, 36] we used 
modified appropriately, not only allow us to perform a semiclassical expansion of 
quantum quantities like the survival probability, as RMT also allows, but give an 
intuitive explanation in terms of classical trajectories. 

In this paper we presented an (almost) complete picture of the semiclassical 
treatment of chaotic systems with tunnel barriers, for arbitrarily complicated 
classically correlated trajectories and for a variety of different quantities. However, 
the examples covered only involved pairs of trajectories, and it is worth remembering 
that when we treat quantities that involve more trajectories additional phases can 
occur for trajectories that never enter the system, as described in [49]. Moreover, 
this picture is only complete, and all the results in this paper are only valid, in the 
particular regime where Ehrenfest time effects are negligible and for times shorter 
than the Heisenberg time. But by placing the possible trajectory structures, their 
contributions and their relations to each other on a firm footing, we hope this work 
could be a stepping stone to a complete expansion based semiclassical treatment of the 
Ehrenfest time regime. Of course the lower order terms have been successfully treated 
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Table Al. The number of trajectory pairs and the number of hnks connecting 



encounters to themselves for syst< 


sms without time reversal symmetry. 


u 






N4(i;) Nb(v) J^6{v) N7(v) 


(2)2 


1 






(3)1 


1 


1 






21 






(2)2(3)1 


49 


5 




(2)1(4)1 


24 




8 


(^) 




4 




(5)1 


8 




8 


(2)« 


1485 






(2)4(3)1 


5445 


189 




(2)-n4)i 


3240 




336 


(2)2(3)2 


4440 


392 




(2)2(5)1 


1728 




420 


(2)1(3)1(4)1 


2952 


168 


392 


(3)3 


464 


84 




(2)1(6)1 


720 




360 


(3)1(5)1 


608 


48 


200 


(4)2 


276 




96 


(7)1 


180 




180 



semiclassically [43, 44, 45, 46], including tunnel barriers [49], and these works show that 
additional trajectory structures appear and contribute in this regime. These additional 
structures and the higher order terms remain to be generalised, but hopefully there is 
some mathematical structure behind these contributions that again can be simplified, 
like in this article, to reproduce continuity and other physical properties. Moving to 
times beyond the Heisenberg time is possibly a much bigger challenge, but progress 
on this front has already been made [24, 25, 26]. 
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Appendix A. Tables of N/(i;) 

Here we record the numbers N;(t;) which are useful for calculating semiclassical 
expansions for transport quantities, like the reflectance, which involve case D. For 
the unitary case, a 2-cncoimtcr cannot return directly to itself, so N2(t') = 0, while 
the remaining numbers are in table Al. The second column was in [18, 19], and is 
only included for reference, while we note that the numbers needed for case C were 
included in [51]. 

For the orthogonal case, we can always add a 2-encounter at the start and end 
of the trajectory so N2(v) = A''(wP~*l), where w^^l is the vector formed from v 
by removing a 2-encountcr. This can be seen, along with the remaining numbers in 
table A2, where the second column is from [18, 19, 51]. The numbers needed for case 
C were included in [51]. 
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Table A2. The number of trajectory pairs and the number of hnks connecting 
encounters to themselves for systems with time reversal symmetry. 



V 


N{v) 


N2(u) 




N4(r) 


N5(v) 


N6(t)) Nt{v) 


(2)1 


1 


1 












5 


1 












4 




4 








(2)3 


41 


5 










(2)1(3)1 


60 


4 


16 








(4)1 


20 






20 






(2)4 


509 


41 










(2)2(3)1 


1092 


60 


132 








(2)1(4)1 


504 


20 




188 






(3)2 


228 




80 








(5)1 


148 








148 




(2)5 


8229 


509 










(2)3(3)1 


23160 


1092 


1592 








(2)2(4)1 


12256 


504 




2388 






(2)1(3)2 


10960 


228 


1968 








(2)1(5)1 


5236 


148 






2392 




(3)1(4)1 


4396 




536 


1164 






(6)1 


1348 










1348 


(2)6 


166377 


8229 










(2)4(3)1 


579876 


23160 


25620 








(2)3(4)1 


331320 


12256 




39448 






(2)2(3)2 


443400 


10960 


48352 








(2)2(5)1 


167544 


5236 






43724 




(2)1(3)1(4)1 


280368 


4396 


19312 


42132 






(3)3 


41792 




8672 








(2)1(6)1 


65808 


1348 








34252 


(3)1(5)1 


52992 




4768 




18016 




(4)2 


24788 






9684 






(7)1 


15104 










15104 



Appendix B. Simplifying channel sum products 

The first sum we wish to simphfy is 

i-i i-i 

fc=l fe=l m 

which just involves geometric progressions. The first part is simply 
/-I 

IJ2pIJ2(^ -Pmf-^ =lJ2Pm [1 - {-^-Pmy-'] , (B.2) 
m fc— 1 m 

while the second is 

/-I 

- i2p'm E - P"^^'~' = - E [1 - (1 - P"^y] + ^ E^^-(i - pr-y~'- 

m k=l m m 

Combining the two parts, the original sum reduces to 
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;-i 

= lG^-Gi=Y,i}-k)h. (B.4) 
fc=i 

We also wish to simplify the following sum 

Z = ki+Y,hGi-k = Y,HkHi-k+i+Y,IkGi-k. (B.5) 

fe=l fe=l /s=l 

The first step is to notice that 

mm m 

= Hk-Ik, (B.6) 
so that we can replace all the Ik terms in (B.5) and obtain 

I i-i i-i 

Z = ^ HkHi_k+i + ^ HkGi-k — ^ Hk+iGi-k 

k=l k=l k=l 

l-l l-l 

= HiHi + ^ Hk {Hi-k+1 + Gi-k) — Hk+iGi-k- (B-7) 

k=l k=l 

The next step is the simplification, 

Hi_k+1 + Gl-k = iPrnil - Pm)'-'' + 1 - (1 - Pm)''-') = ^ (l - (1 - PmY''') 
m m 
= Gi-k+1, (B-8) 

which we can put into (B.7). We also change the sum index on the second sum (with 
k' = k + 1), to obtain 

;-i / 
Z = HiHi + Y^HkGi-k+i - ^ Hk'Gi-k'+i 

k=l fe'=2 

= HiHi + HiGi - HiGi, (B.9) 

since all the terms in the two sums mutually cancel apart from the start and end ones. 
As Hi = Gi, we obtain the final result of 

/-I 

Z = GiGi = ki + Y,IkGi-k. (B.IO) 
fc=i 
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